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I.  Introduction. 


The  research  described  below  was  completed  during  the  period  January  1,  1988,  to 
December  31,  1990.  This  work  was  funded  by  the  Air  Force  of  Scientific  Research  under 
grant  AFOSR-SS-0074. 

In  recent  years  the  feasibility  and  advantages  of  active  control  surfaces  to  reduce  ma¬ 
neuver,  gust,  and  fatigue  loads,  and  dampen  vibration  that  contributes  to  flutter  have  been 
extensively  studied.  A  systematic  procedure  for  control  design  requires  the  development  of 
a  realistic  mathematical  model  that  predicts  the  dynamical  behavior  of  the  physical  system. 
The  development  of  state  space  models  for  aeroelastic  systems,  including  unsteady  aerody¬ 
namics,  is  potentially  import  for  the  design  and  development  of  highly  maneuverable  aircraft. 
A  goal  of  this  research  effort  was  to  develop  a  complete  dynamical  model  for  an  aeroelastic 
system  in  which  the  elastic  motions  of  the  structure  are  coupled  with  the  motion  of  the 
surrounding  fluid.  Our  study  included  the  task  of  identifying  the  appropriate  state  spaces 
for  which  the  corresponding  mathematical  model  was  well-posed  as  well  as  the  development 
and  analysis  of  computational  algorithms  for  these  systems. 

It  has  been  shown  that  fractional  order  operators  can  provide  accurate  models  for  vis- 
coelasticiallv  damped  structures.  The  governing  equation  for  such  systems  has  the  form 
of  a  partial-integral  equation.  These  systems  can  be  viewed  as  neutral  functional  differ¬ 
ential  equations  having  a  structure  common  to  the  mathematical  models  described  above 
for  aeroelastic  systems.  For  our  study  of  viscoelastic  systems  we  focused  our  attention  to 
the  task  of  developing  a  parameter  identification  scheme  for  these  systems.  As  a  first  step, 
we  investigated  the  effectiveness  of  quasilinearization  for  parameter  identification  in  both 
linear  and  nonlinear  parabolic  partial  differential  equations.  Our  study  included  the  heat 
equation  where  the  unknown  parameter  was  the  spatially  varying  diffusion  coefficient  and 
the  nonlinear  Burgers’  equation. 
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II.  Aeroelastic  Systems. 

A  completed  dynamic  model  for  the  elastic  motions  of  a  three-degree-of-freeduin  “typ¬ 
ical”  airfoil  section,  with  flap,  in  a  two-dimensional  incompressible  flow  was  formulated  by 
Burns,  Cliff  and  Herdman,  (see  [1,  reference  S]),  in  terms  of  a  functional  differential  equation 
of  neutral-type.  In  subsequent  papers,  Burns,  Herdman  and  Turi  (see  [1,  references  10-13 
and  24]),  the  well-posedness  of  the  modeling  neutral  equation  was  studied  in  a  product 
space  framework.  The  analysis  showed  that  the  dynamic  model  extends  to  a  well-posed 
state-space  model  on  the  product  space  IR8  x  Lp\  p  >  1  if  and  only  if  p  <  2.  Since  the 
ultimate  goal  of  the  modeling  process  was  to  generate  a  framework  for  the  design  of  active 
control  schemes  for  flutter  suppression  these  results  would  suggest  that  one  consider  histories 
belonging  to  Lpi  p  <  2.  On  the  other  hand  in  the  derivation  of  the  evolution  equation  for 
the  circulation  on  the  airfoil  one  has  to  assume  that  the  circulation  history  belongs  to  Lp  for 
p  >  2  in  order  to  guarantee  the  applicability  of  Sohngen’s  inversion  formula  for  obtaining 
a  representation  for  the  solution  of  the  airfoil  equation.  Motivated  by  these  observations, 
we  investigated  the  problem  of  finding  a  state-space  such  that  Sohngen’s  inversion  formula 
was  applicable  and  at  the  same  time  the  resulting  neutral  system  was  well-posed.  That  is. 
we  wanted  to  identify  conditions  on  the  initial  data,  including  the  past  history,  to  guaran¬ 
tee  well-posedness  of  the  model  equations  that  matched  the  conditions  necessary  to  justify 
the  validity  of  the  inversion  formula  employed  to  obtain  the  solution  of  the  airfoil  equation. 
During  this  time  period  we  have  shown  that  this  compatibility  can  be  obtained  by  assuming 
that  the  circulation  history  belongs  to  a  weighted  L2  space. 

Our  study  of  the  well-posedness  of  functional  differential  equations,  including  the  aeroe¬ 
lastic  system  discussed  above,  has  been  in  the  context  of  functional  analytic  semigroup 
theory.  The  applicability  of  the  semigroup  approach  for  the  development  of  approximation 
techniques  for  parameter  estimation  and  optimal  control  for  such  system  requires  the  iden¬ 
tification  of  an  appropriate  state-space  which  not  only  yields  the  Co-semigroun  hut  Hso 
allows  one  to  obtain  a  dissipative  estimate  for  the  infinitesimal  generator  for  the  semigroup. 


This  requirement  is  due  to  the  convergence  of  the  numerical  schemes  being  established  by 
using  the  Trotter-Kato  semigroup  approximation  theorem.  This  dissipative  estimate  for  the 
aeroelastic  system  with  state-space  IRS  x  L?  is  guaranteed,  however,  finding  the  equivalent 
norm  on  that  space,  which  is  needed  to  give  the  dissipative  estimate,  is  not  straightforward. 
For  retarded  functional  differential  equations  and  atomic  neutral  functional  differential  equa¬ 
tions  on  the  state  space  IR"  x  Li,3 ,  g  a  weight  function,  a  well-developed  theory  is  available. 
However,  for  singular  neutral  systems,  such  as  our  aeroelastic  system,  a  characterization  of 
g  that  is  needed  to  obtain  the  dissipative  estimate  was  unknown.  Our  investigation  shows 
that  the  selection  of  g  plays  a  key  role  in  the  development  of  the  semigroup  formulation  for 
our  system. 

In  [2]  we  establish  that  the  finite  delay  version  of  the  aeroelastic  system  generates  a 
Co-semigroup  on  the  product  space  IR'  x  £2,3  and  show  that  the  infinitesimal  generator  of 
that  semigroup  does  indeed  satisfy  a  dissipative  estimate.  Here  the  function  g  is  defined 
by  g($)  =  [1  +  (— s)"1]2.  Our  framework  can  be  extended  to  include  more  general  singular 
neutral  functional  differential  equations.  In  the  general  case  we  allow  a  weight  function  g 
satisfying  (i)  g(s)  >  0,  g(-)  €  Ti,  g  >  0  and  (ii)  for  any  given  c  >  0,  there  is  a  r  >  0  such 
that  g(s)  >  cg(s)  on  [— r,  0].  In  [3]  we  present  an  approximation  scheme  and  provide  two 
illustrative  examples  to  indicate  the  feasibility  of  our  approach.  One  example  is  a  scalar 
equation  with  inconsistent  initial  data,  no  smoothness  on  the  solution  at  t  =  0,  the  other 
example  is  a  two-dimensional  system  which  exhibits  characteristics  similar  to  those  one 
observes  in  the  aeroelastic  system  (see  [2,3]). 

In  [1]  we  study  the  derivation  of  the  equations  for  the  aeroelastic  system  assuming  that 
the  circulation  history  belongs  to  a  weighted  £,2  space.  We  establish  that  the  assumption  that 
the  “past  history”  of  the  derivation  of  the  total  circulation  function  belongs  to  Li,g,  same 
g  that  provided  the  dissipative  estimate,  is  sufficient  to  assure  to  validity  of  the  Sohngen's 
inversion  formula  for  the  solution  of  the  airfoil  equation.  Therefore,  we  have  obtained  com¬ 
patibility  of  the  validity  of  the  inversion  formula  and  the  well-posedness  of  the  resulting 
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neutral  system.  The  resulting  state  space  IR'  x  Li,3  and  the  numerical  techniques  described 
above  provides  a  suitable  setting  for  control  design  for  the  aeroelastic  system. 

For  the  complete  system  described  above  the  evolution  equation  for  the  circulation  on 
the  airfoil  was  coupled  to  the  rigid-body  dynamics  of  the  airfoil  to  obtain  the  infinite  delay 
singular  neutral  system  of  the  form 

d.  f° 

—  Ax(t)  +  /  /l(s)  x  (t  4  s)  ds  =  Bx(t)  4  /  B(s)  x  (t  4  s)  ds  +  f(t)  (2.1) 

«  L  2-00  2-00 


as  the  mathematical  model.  In  (2.1)  the  S  x  S  matrix  A  is  singular,  the  8  x  S  matrix  function 
A(s)  is  weakly  singular,  <488(,s)  =  {{Us  —2)/Us)*  and  /  represents  a  forcing  term  (possible 
control)  for  the  system.  The  state  of  the  system  (2.1)  includes  the  past  history  of  T  where 
T  denotes  the  total  airfoil  circulation  (see  [1]).  This  past  history  of  T  may  not  be  observable 
for  the  entire  past  time  (— oo,0).  However,  one  would  be  able  to  observe  this  past  history- 
over  a  finite  time  interval  say  [— r, 0].  Since  the  kernel  function  A(s)  is  not  integrable  on 
(—oo,0)  it  is  not  possible  to  address  the  infinite  delay  problem  as  a  finite  delay  problem  by 
mapping  the  interval  (— oo,0)  to  some  finite  interval  [— T,  Oj.  On  the  other  hand  the  ^finite 
delay”  version  of  (2.1) 

=  Bx{t)+  f  B{s)  x  (t  +  s)ds  +  f{t)  (2.2) 


d_ 

dt 


Ax{t)  4- 


/> 


s)  x  {t  +  s)  ds 


as  mention  earlier  has  been  extensively  studied.  During  this  period  we  have  investigated  the 
possibility  of  taking  advantage  of  our  results  for  the  finite  delay  system  to  study-  the  infinite 
delay  system.  In  [4]  we  develop  an  approximation  technique  for  (2.1)  on  an  interval  of  the 
form  [0,T]  using  the  approximation  techniques  we  developed  for  the  finite  delay  system 
(2.2).  In  order  to  accomplish  this  we  have  made  use  of  the  special  structure  found  in  the 
aeroelastic  system.  In  particular,  the  representation  of  the  aeroelastic  system  allows  one  to 
view  the  system  as  having  an  integrodifferential  Volterra  component  and  a  singular  integral 
component  (see  [4]).  We  have  shown  that  the  infinite  delay  singular  integral  component  found 
in  (2.1)  can  be  approximated  using  the  corresponding  finite  delay  singular  integral  component 
of  (2.2).  Our  results  including  error  estimates  for  the  approximation  aie  presented  in  [4], 
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Our  two  step  approximation  procedure  starts  with  truncating  the  infinite  delay  term.  The 
approximation  is  dependent  on  the  specified  initial  data.  The  desired  error  bound  together 
with  the  time  interval  [0,  T]  on  which  one  wishes  to  find  the  solution  dictates  the  necessary 
truncation.  The  second  step  consists  of  employing  the  approximation  scheme  established  for 
the  associate  finite  delay  system. 


III.  Parameter  Identification. 

During  this  period  we  developed  a  scheme  for  identification  of  unknown  parameters  in 
parabolic  partial  differential  equations.  The  task  that  we  undertook  was  to  develop  a  quasi¬ 
linearization  based  algorithm  that  could  be  used  to  identify  parameter  for  nonlinear  as  well 
as  linear  parabolic  partial  differential  equations.  As  a  first  step  we  considered  a  parameter 
identification  for  the  heat  equation  where  the  unknown  parameter  was  the  spatially  varying 
diffusion  coefficient.  Our  study  included  theoretical  detail  and  careful  numerical  testing. 
V/c  established  the  Frechet  differentiability  of  the  solution  with  respect  to  the  unknown  pa¬ 
rameter  which  was  shown  to  be  a  primary  result  necessary  to  ensure  local  convergence  of 
the  approximation  scheme  used  to  find  the  unknown  parameter.  The  numerical  examples 
that  we  present  demonstrate  the  rapid  convergence  and  accuracy  of  the  algorithm.  Our  ef¬ 
forts  to  date  involving  the  use  of  quasilinearization  in  nonlinear  partial  differential  equations 
have  only  been  numerical  ones.  However,  the  success  of  our  numerical  testing  leads  us  to 
believe  that  our  quasilinearization  scheme  will  provide  an  effective  approach  to  parameter 
identification  for  certain  nonlinear  equations.  Our  results  are  presented  in  [5], 
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IV.  Interaction  with  Air  Force  Laboratories  and  Industry. 

During  the  peiiod  of  this  grant,  Dr.  Herdman  has  consulted  with  Air  Force  Laboratory 
and  Industry  personnel  for  motivation  and  guidance  for  this  research  effort.  We  discussed 
our  research  results  and  possible  implementation  with  the  following  people: 

•  Wen-Huei  Jou,  Manager,  CFD  Development,  Computational  Fluid  Dynamics  Lab.,  Boe¬ 
ing  Commercial  Airplanes,  Seattle,  Washington. 

•  Lewis  F.  Jurey,  Manager,  Flutter  and  Dynamics,  Advanced  Development  Projects,  Lock¬ 
heed  Aeronautical  Systems  Company,  Burbank,  California. 

•  J.  Jenkins  and  J.  Lee,  Wright  Research  and  Development  Center,  Wright-Patterson 
AFB,  Ohio. 

V.  Personnel  Supported. 

Principal  Investigator:  Terry  L.  Herdman-Professor,  Department  of  Mathematics. 
Graduate  Project  Assistant:  Patricia  W.  Hammer-Ph.D.  Student,  Department  of  Mathe¬ 


matics. 
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submitted. 
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